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ABSTRACT 



Using data from the DEEP2 Galaxy Redshift Survey and HST/ ACS imag- 
ing in the Extended Groth Strip, we select nearly 100 interacting galaxy systems 
including kinematic close pairs and morphologically identified merging galaxies. 
Spitzer MIPS 24 yum fluxes of these systems reflect the current dusty star for- 
mation activity, and at a fixed stellar mass (M^, ) the median infrared luminosity 
{LiR ) among merging galaxies and close pairs of blue galaxies is twice (1.9±0.4) 
that of control pairs drawn from isolated blue galaxies. Enhancement declines 
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with galaxy separation, being strongest in close pairs and mergers and weaker in 
wide pairs compared to the control sample. At J ~ 0.9, 7.1% ± 4.3% of massive 
interacting galaxies (M* > 2 x lO^^M©) are found to be ULIRGs, compared to 
2.6% ± 0.7% in the control sample. The large spread of Ljji/M^ among inter- 
acting galaxies suggests that this enhancement may depend on the merger stage 
as well as other as yet unidentified factors (e.g., galaxy structure, mass ratio, 
orbital characteristics, presence of AGN or bar). The contribution of interacting 
systems to the total IR luminosity density is moderate (< 36%). 

Subject headings: galaxies:interactions - galaxies:evolution - large-scale structure 
of universe - infrared: galaxies 



1. INTRODUCTION 



Galaxy-galaxy interaction has long been regarded as a key process in galaxy evolution, 
espec i ally as a mechanism for enhancing star formation during mergers (ILarson fc Tinsley 



19781 : iBarton et al.ll2000l ). Hydrodynamic A^-body simulations show that active star for- 
mation can be triggered by gaseous i i iflows resu lting from mergers of gas-rich galaxies 
( iMihos &: Hernquistlll996l : iBarnej |2004J : ICoxll2004l ). Interaction-triggered star formation is 
also thought to be responsible for luminous infrared sources. In the local universe, luminous 
infrared gala xies (LIRGs) and ultra lurninous infrared galaxies (ULIRGs) are primarily merg- 
ing systems (ISanderslll988l : iBorne et al.lll999l ). Nevertheless, by studying various star forma- 
tion in dicators of interacting galaxies and normal galaxies, iBergvall. Laurikainen. fc Aalto 
( 120031 ) concluded that galaxy interactions in general are inefficient triggers of starbursts; 
interactions are a necessary but not sufficient condition to trigger violent starbursts. 

The importance of galaxy interactions in the volume-averaged galaxy star formation rate 
(SFR) remains an open question. Recent studies of mid-IR (MIR) sources at a median red- 
shift of z ~0.7 suggest that the IR density at that ep och is dominated by morphologically nor- 

mal g alaxies instead of strongly interacting galaxies ( Bell et aPboOS : Melbourne. Koo. &: Le Floc'h 
20051 ). Two main factors may contribute to this result. First, only a small fraction of the 



galaxy population may b e undergoing a major merger at any given time ( ICarlberg et al. 



Bundv et al. 


2004; 


Lin et al. 


2004) 



z ~0.7 may be enhanced relative to the local population (iBell et al.ll2005l ). perhaps as the 
result of internal processes such as a higher gas fraction leading to a higher SFR, such 
that galaxy interactions may have less dramatic effects on triggering star formation at that 
epoch and/or may be harder to identify. It is the aim of this Letter to examine this second 
hypothesis. 
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This Letter presents an analysis of the IR properties of close kinematic galaxy pairs, 
morphologically selected merging galaxies, and a control sample of randomly selected pairs of 
isolated galaxies, in the range 0.1 < z < 1.1. In §3 we show our analysis of the IR luminosity 
[LiR ) versus stellar mass (M^, ) for the interacting galaxies and control samples, and the 
relation between the IR luminosity-to-mass ratio {Lir/M^ ) and the projected separation 
of the galaxy pairs. Discussion and conclusions are given in §4. Throughout this Letter we 
adopt the following cosmology: Hq = lOO/i km s~^Mpc~^, VL^ = 0.3, and = 0.7; h = 0.7 
is adopted when calculating the rest-frame magnitude. Magnitudes are given in the Vega 
system. 



2. DATA, SAMPLE SELECTIONS, AND METHODS 



2.1. Data 



The sample used in this Letter consists of ~ 4000 optical selected galaxies with secure 
redshifts from the DEEP2 Galaxy Redshift Survey and MIPS imaging in the Extended Groth 
Strip (EGS) region. The EGS has extensive multi-wavelength observations from both the 
ground and spa ce, coordinated by the All- Wavelength Extended Groth International Survey 
(AEGIS) team JPavis et al.ll2007[ ). The spectral resolution of the DEEP2 survey is R ~ 5000, 
corresponding to redshift errors of ~30 km s~^, allowing us to define kinematic pairs (see 
§ 2.2). The average sampli ng rate of the slit pla ced on galaxies is ~ 60%, and the average 
redshift success is ~ 73% (IWillmer et al.ll2006l ). Deta ils of the survey a nd /T-correction 



procedure are described in a series of DEEP2 papers (IDavis et al.l l2003l : I Coil et al.l 12004 



Willmer et al.ll2006l ). The MIPS (Multiband Imaging Photometer for Spitzer; Rieke et al. 



2004) 24 /im observations used here were carried out on 2004 June 19 and 20; iDavis et al. 
( 120071 ) provide further details. 



2.2. Selection of Kinematic Pairs, Merging Galaxies, and Control Samples 

Three classes of galaxy samples are investigated in this study: kinematic pairs, merg- 
ing galaxies, and a control sample of randomly selected pairs of isolated galaxies, in three 
redshift bins in the range 0.1 < z < 1.1. To select kinematic galaxy pairs, we use the 
following criteria: (1) For each galaxy in the spectroscopic redshift sample, we first search 
for any kinematic companions with a relative line-of-sight velocity Av < 500 km s~^and a 
projected physical separation (onto the plane of the sky) Ar^ <300 h'^kpc . (2) Among these 



companions, close pairs are identified such that they satisfy AVp < 50 h ^kpc (jPatton et al. 
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20021 : iLin et al.ll2004l ). We find a total of 56 close pairs. (3) For comparison, wide pairs are 
identified as kinematic companions with 50 h'^pc < AVp < 300 h~^kpc , and we retain 
only wide pairs that are closest companions to each other (such that there is only one wide 
pair companion per galaxy). The spectroscopic redshift sample is not entirely complete, 
however; it is possible that for a given galaxy, there exists a closer companion that is missing 
in the redshift sample. Therefore, AVp for each paired system is an upper limit on the actual 
distance to the closest companion. To minimize this effect we also search for photometric 
companions of each galaxy, and we keep only those wide pairs that do not have a photometric 
companion (with observed Am/j < 1 mag) within AVp < 30 /i^^kpc , for a sample of 126 wide 
pairs. (4) Absolute magnitude limits were applied to all pair samples within a given redshift 
bin (see Table 1). (5) To study the effect of interactions between two gas-rich galaxies, we 
further applied a color cut using rest-frame U — B = 0.25 as the division between blue and 
red galaxies and only include blue galaxies in our samples. It is possible that we have missed 
some red, dusty, late-type galaxies due to our color selection. Howe ver, the number de nsity 



of that population is ~8% compared to the blue galaxy population (IWeiner et al.ll2005l ) and 
is therefore negligible for our analysis. 

We also define a set of morphologically identified merging galaxies using a subsample 
of blue galaxies in the EGS that have deep HST images taken with Advanced Camera 
for Surveys (ACS) as part of GO program 10134 (PI: M. Davis). These morphologically 
identified merging galaxies are preferentially in strongly interacting systems and are in a 
different merger stage than the kinematic close pairs. We identify merging galaxies using 
three non-parametric parameters: the Gini coefficient (G), the second-order moment of 
the brightest 20% of a galaxy's pixels (M20), and th e asymmetry measurement (A), all of 



which have been shown to ident ify merger candidates (IConselice II2003I : IConselice et al.ll2003 



Lotz. Primack. fc Madaull2004j ). We first select galaxies with morphological parameters of 
G > -0.115 * M20 + 0.384 OT A > 0.25 (see Lotz et al. 2007 for discussion on G and 
M20 in detail), and then we perform by-eye examinations of each object to keep only those 
candidates that show apparent interaction signatures (eg., tidal tails, distorted morphology, 
and double nuclei). This results in a sample of 56 merging systems. 

For a fair comparison with isolated galaxies we also construct a sample of 1800 control 
pairs, each of which consists of two galaxies randomly selected from galaxies that are iso- 
lated, with no spectroscopic companion within 100 /i^^kpc or photometric companion (with 
observed Am^j < 1 mag) within 30 /i^^kpc . Isolated galaxies have the same magnitude and 
color cut as adopted for interacting galaxies. 
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2.3. Stellar Mass and Total IR Luminosity 



Stellar masses are deri ved from rest-fra r ne (B — V) colors and absolu te Mb magnitude as 
described by the models of iBell &: de Jond (I2OOII ) and iBell et al.l ( 120051 ) . who find a scatter 
of ~0.3 dex in the resulting stellar mass estimates. These measuremeri ts were further refined 
by comparison with the stellar masses derived by iBundy et al.l (120061 ) using detailed fits to 
the spectral energy distribution (SED). Empirically we find that the difference between the 
masses estimated using the rest-frame colors and the SED fits is improved by making small 
corrections for the redshift, rest-frame {U — B) and {B — V) colors: 



log 



10 



Me 



-OA{Mb - 5.48) + 1.737(5 - V) + 0.098(f/ - B) 
-0.130(t/ - Bf - 0.2682 - 1.003. 



These terms have the effect of correcting the z = measurements of iBell fc de JongI (120011 ) 
to the galaxy redshift as well as accounting for evolution in color. Masses estimated from 
Equation 1 agree with those from full SED fits (when available; Bundy et al. 2006) to an 
rms accuracy of 0.25 dex. 

Local studies have shown that the rest-frame MIR luminosity is tightly correlated with 
the to tal IR luminosity over a wide range of galaxy types ( Roussel et al.lboOl : Papovich &: Bell 



20021). This relation appears to hold up reasonably well to z ~ 1 (lElbaz et al 



2002 



Appleton et al.l |2004| ) and hence allows us to estimate the IR luminosity from the ob- 
served flux of 24 fim corresponding to the rest-frame 11 -21 nm flux over the re dshift range 
0.1 < z < 1.1. Following the same procedure adopted by lLe Floc'h et al.l (120051 ) . we convert 
the o bserved 24 nm fl ux in to the total IR luminosity by fltting for each source using SEDs 
from ICharv fc ElbazI J200lh . Q Because of the large PSF (FWHM 6") of the 24 /xm data, 
we work on the total IR luminosity in pairs rather than the IR luminosity of each individual 
galaxy. This also allows us to ma ke comparisons to local ULIRGs , which are found to be 
highly distorted merg ing systems (ISanderslll988l: Borne et al.lll999l ). and to merger simula- 
tions in the literature (lBarnesll2004j : I Jonsson et al.ll2006l ). which yield the total star formation 
rate of the merger system. For each pair, we then calculate the mean IR luminosity-to-mass 
ratio, LiFt/M^, , by summing the IR luminosity of both galaxies and dividing by the total 
stellar mass of both components. 



Various template libraries lead to a scatter in this conversion of a factor of 2-3 |Le Floc'h et al-lboosh . 
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3. RESULTS 

Figure 1 shows the relation between Lm and M^^ for three of the galaxy samples. A 
relatively tight correlation is seen for the control pairs, and the close pairs and mergers 
are within the bounds of this relation but tend to occupy the upper region of L/^ for a 
given Mj, . In the highest redshift bin, the fraction of merging galaxies and close pairs more 
massive than 2 x IO^^Mq that are ULIRGs is 7.1% ± 4.3%, compared to 2.6% ± 0.7% for 
the control sample. To address quantitatively the difference of Ljr among those samples, we 
compute the median Lm/M^ for each sample. Using the median value of L/ij/M^, rather than 
performing the K-S test avoids the problem of handling the upper limits. When calculating 
the median Ljji/M^ , stellar mass cuts are further applied in each redshift bin (see Table 1). 
Over the entire redshift range, the median L/ij/M* of pairs and mergers is twice (1.9 ± 0.4) 
that of the control sample. To assess if our results are sensitive to the method we have 
adopted to estimate stellar masses, we also apply the same analysis on a subs ample of the 



data w ith stellar mass measurements derived from the SED-fitting procedure of iBundy et al. 



( 120061 ). The enhancement in L/^j/M^, increases to 3.1 ±1.3 in this case, consistent with our 



initial result within the la uncertainties. 

An alternative way to investigate the effect of galaxy interactions on the star formation 
rate is to study the dependence of Lj^i/M^, on the pair separation, as shown in Figure 2. For 
comparison, distributions of Lm/M^, in control pairs are also shown along the right-hand 
axes. The merging galaxies and close pairs possess higher median Lir/M^, than wide pairs 
and control pairs, as shown in Table 1, but also have a wider spread in Lm/M^ . The 
declining envelope of Ljj^/M^, as a function of Ar^ is similar to the behavior of the local pair 



results reported by iBarton et al.l (120001 ) . who use various sets of emission lines as the star 
formation tracer. The median L/^/M* of wide pairs, on the other hand, becomes close to 
that of the control pairs, indicating that the effect of galaxy interaction on star formation 
activity is limited to several tens of /i^^kpc . 



4. DISCUSSION AND CONCLUSION 



Using data from the DEEP2 Galaxy Redshift Survey with MIPS 24 /xm imaging in the 
EOS, we find that the combined IR luminosity at a given stellar mass of blue merging galaxies 
and kinematic close pairs is greater by a factor of 1.9 ± 0.4 than that of control pairs ran- 
domly drawn from blue isolated ga laxies. This enhancement is consistent with low-redshif t 
studies of the SFR in galaxy pairs (ILambas et al.ll2003l : iNikolic. CuUen. fc Alexander l2004j ). 
In addition, we also observe a declining envelope of the Ljji/M^ with increased projected 
separation of kinematic galaxy pairs. Based on the assumption that the IR emission is tightly 
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associated with the SFR, our resuhs quahtatively support the picture of tidally triggered 
starbursts as predicted in hydrodynamic simulations, ahhough the effect is apparently not 
large. The frequency of ULIRGs in massive interacting systems is only 7.1% ± 4.3% in the 
highest redshift bin (0.75 < z < 1.1), and no ULIRGs are found below a combined stellar 
mass for the pair of 4 x 10^^ Mq. The mass ratio between two g as-rich galaxies may be a key 
element in generating ULIRG luminosities (jPasyra et al.ll2006l ). along with the dependence 
on the combined stellar mass of the interacting system. The wide spread of specific star 
formation rates in interact ing galaxies found h ere (as inferred from the broad distribution of 
Lir/M^, ), and also seen in lBarton et al.l (|2000| ). indicates that some other mechanism(s) may 
determine the strength of the induced star formation as well - e.g., the galaxy structure, the 
mass/luminosity rat io, the orbital config urati on, and the existence of AGN aii d/or bar. For 
example, studies by iLambas et al.l (120031 ) and lWoods. Geller. fc BartonI (120061 ) suggest that 
pairs with comparable luminosity do show stronger star formation activity than pairs with a 
larger luminosity contrast. The scatter in Lir/M^, that we find could also be due to the fact 
that these systems are in different stages of merging: some kinematic pairs are likely on their 
first approach while others are being seen after their first passage. Another possibility is that 
some of the close pairs selected by our Ar^ and Af criteria are not physical close pairs in 
real space (as opposed to redshift space; Perez et al. 2006). Our results nevertheless provide 
a constraint on the amount of induced star formation activity due to galaxy interactions at 
z ~ 0.1 - 1.1. 

Finally, we note that tidally triggered star formation contributes moderately to the 



high IR luminosity density at intermediate redshifts of 0.4 < z < 1.1. iLin et al.l (120041 ) 
estimated that the pair fraction (Arp <50 /i~^kpc and Av <500 km s~^) of L* galaxies over 
this redshift range is less than 15% and that the redshift evolution in the pair fraction is also 
much lower than the rapid decline of IR luminosity density seen with decreasing redshift. 
Additionally, the fraction of r norphologically id entified merging galaxies remains roughly 
constant at 7% up to 2 ~ 1 (ILotz et al.l l2006l ). These results, when combined with the 
moderate increase of star formation activity seen in interacting systems shown here, indicate 
that the contribution of galaxy interactions to the total IR density at intermediate redshift is 
< 36%. This is consis tent with the con clusi ons from the studies of IR population s of various 
morphology types by iBell et al.l (120051 ) and iMelbourne. Koo. fc Le Floc'hl (120051 ) . 
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Table 1. Properties of the pair sample 



z 


Sample 


Mb (individual) 


M* (combined) 


No. 


^24 


median Lzr/M^ 


0.1 < z <0.4 


MG 


-21.7 < Mb < -18.7 


9.7<logioM,<10.7 


3(9) 








close pair 


-21 < Mb < -18 


9.7<logioM.<10.7 


10(10) 


68±10(%) 


2.2±0.7 




wide pair 


-21 < Mb < -18 


9.7<logioM*<10.7 


17(20) 


89±6(%) 


1.3±0.1 




control pair 


-21 < Mb < -18 


9.7<logioM*<10.7 


480(600) 


74±11(%) 


1.0±0.3 


0.4 < 2 <0.75 


MG 


-22.7 < Mb < -19.7 


10<logioM,<ll 


13(18) 


62±13(%) 


3.1±0.8 




close pair 


-22 < Mb < -19 


10<logioM.<ll 


21(24) 


71±7(%) 


2.7±0.4 




wide pair 


-22 < Mb < -19 


10<logioM*<ll 


44(54) 


64±5(%) 


1.7±0.3 




control pair 


-22 < Mb < -19 


10<logioM*<ll 


502(600) 


60±7(%) 


1.5±0.4 


0.75 < 2 <1.1 


MG 


-23.7 < Mb < -20.7 


10.3<logioM*<11.3 


20(29) 


73±10(%) 


3.3±2.7 




close pair 


-23 < Mb < -20 


10.3<logioM,<11.3 


22(22) 


76±7(%) 


3.8±0.8 




wide pair 


-23 < Mb < -20 


10.3<logioA^,<11.3 


45(52) 


62±5(%) 


1.5±0.2 




control pair 


-23 < Mb < -20 


10.3<logioM.<11.3 


547(600) 


64±7(%) 


1.9±0.5 



Note. — Mb is the rest-frame B-band magnitude of each galaxy in pairs. The luminosity range of merging galaxies (MGs) 
is chosen to be twice as bright as that of individual galaxy in pairs because merging galaxies are assumed to be comprised 
of two galaxies. M* is the stellar mass cut when computing the median L/^j/M* and 24 jum detection rate d24- In the fifth 
column, the values outside (inside) the parentheses axe the numbers of paired systems with (without) a stellar mass cut. 
Error bars quoted here are computed by bootstrapping. 
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Fig. 1. — Total IR luminosity (L/_r ) vs. stellar mass (M* ) for kinematic close pairs 
(blue circles), control pairs (green crosses), and morphologically identified merging galaxies 
(magenta triangles). Systems without a 24 /xm detection are assigned low Ljr values < 3x 10^ 
Lq. Horizontal red solid lines represent the IR thresholds for ULIRGs [Ljr > 10^^ Lq) and 
LIRGs {Ljji > 10^^ Lq). Black dotted lines correspond to 70 /iJy (Sex) in 24 fj,m at lower 
and upper redshift limits. Arrows in the lower right corner of each panel show the scatter 
(0.3 dex) in the conversion of 24 /im flux into total between different template libraries. 
It can been seen that the close pairs and merging galaxies possess higher Lm than control 
pairs for a given . 
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Fig. 2. — Lir/M^ as a function of projected separation for kinematic galaxy pairs (blue 
circles) and merging galaxies (magenta triangles) in three redsliift bins. Merging galaxies 
are assigned to Ar^ ~ 5 /i~^kpc , which is the minimum separation of kinematic close pairs. 
Galaxies with Lm > 10^^ Lq (ULIRGs) are shown as larger symbols. Distributions of 
Ljji/M^: in control pairs are also shown along the right-hand axes for comparison. Data 
points at Ljji/M^ = refer to those sources with no 24 /im detection. The number within 
the parentheses in the bottom panel denotes the L/ij/M* of the closest data point. Close 
pairs (Arp < 50 h~^kpc ) and merging galaxies are found to have a higher median Lm/M^ 
and a wider spread of Li^/M^ than those of wide pairs (Ar^ > 50 /i^^kpc ) or control pairs 
(see Table 1). 



